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ABSTRACT 

We perform a three flavour analysis of the atmospheric, accelerator and reac- 
tor neutrino data from the Kamiokande, LSND and Bugey experiments respectively. 
Choosing the values of Am 2 obtained from two flavour fits of the first two experiments, 
the allowed ranges of the three generation mixing angles are determined. The acceler- 
ator experiments CHORUS and NOMAD are found to be most sensitive to regions of 
the allowed parameter space which correspond to genuine three generation solutions 
for the atmospheric neutrino anomaly. 
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In the standard model of electroweak theory the neutrinos are assumed to be 
massless. But there is no compelling theoretical reason behind this assumption. If 
the neutrinos are massive then, as in the quark sector, the weak interaction basis of 
neutrinos may be different from the mass eigenstate basis - leading to mixing between 
different flavours. A way for probing such mixing and small neutrino masses is provided 
by neutrino oscillations. Two well known neutrino puzzles that can be explained by 
flavour oscillation of neutrinos are the solar neutrino problem and the atmospheric 
neutrino anomaly. The recent declaration by the Liquid Scintillator Neutrino Detector 
(LSND) collaboration |l| that they are observing an excess of I7 e s (over the expected 
background) which can be attributed to V^—V e oscillations has added a new dimension 
to the issue of neutrino mass and mixing. LSND is most sensitive to Am 2 ~ 6eV 2 
and the significance of this result for particle physics, astrophysics and cosmology 
has been investigated || f|]. One notes that the three phenomena mentioned above - 
namely, the solar neutrino problem, the atmospheric neutrino anomaly and the — 
V e oscillations observed by the LSND group - require vastly different mass ranges. 
The solar neutrino problem can be explained either by Mikheyev-Smirnov-Wolfenstein 
oscillation B for Am 2 ~ 6 x 10 _6 eV A2 and sin 2 29 ~ 7 x 10~ 3 (non-adiabatic solution) 
and Am 2 ~ 9 x !Q~ 6 eV 2 and sin 2 29 ~ 0.6 (large mixing angle solution) || or by 
oscillation in vacuum for Am 2 ~ (0.45 - 1.2) x lO -10 ^ 2 and sin 2 29 ~ (0.6 -1.0) J7j in 
a two generation scenario. The atmospheric anomaly can be explained by either v ^ — u e 
or Ufj, — v T oscillations in a two generation picture. The analysis of the new multi-GeV 
data as well as the previous sub-GeV data of the Kamiokande collaboration predicts the 
following best-fit parameters (Am 2 , sin 2 29) = (1.8 x lO~ 2 eV 2 , 1.0) for v li —v e oscillation 
and (1.6 xlO^eV 2 , 1.0) for z/ M - v T oscillation @. 

Although each experiment can be explained by two flavour neutrino oscillations, 
there are several motivations to go beyond this approximation. The LEP result that 
there are three light neutrinos is also supported by the requirements of nucleosynthesis 
in the early universe. In the quark sector, mixing between three generations is well 
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established. A natural question then is how do experiments constrain three neutrino 
mixing? We stress that in a realistic three flavour framework it is important to do a 
combined analysis to find out the allowed range of parameters rather than using sep- 
arate two flavour schemes. In particular, this might uncover regions in the parameter 
space sensitive to the presence of the third generation which cannot be probed in the 
two flavour limit. 

In this paper we perform a three flavour analysis of the atmospheric and LSND 
data assuming that the presently reported values will not change significantly as more 
results accumulate. The constraints obtained from the reactor experiments are also 
incorporated. We take the Am 2 s as: A 12 ~ A 13 = 6eV 2 in the LSND range and 
A 2 3 = lQ~ 2 eV 2 as preferred by the atmospheric neutrino data. The ~ sign means 
we neglect 10~ 2 as compared to 6. It will become clear as we proceed that most of 
our analysis does not depend on this specific choice as long as the order of magnitude 
remains the same. The three mixing angles are allowed to cover the whole range from 
to 7r/2. For atmospheric neutrinos, we find, in addition to the two flavour results, 
genuine three generation solutions where both — v e and — v T oscillation channels 
simultaneously contribute. The implications of the allowed areas thus obtained for the 
accelerator experiments CHORUS and NOMAD searching for z/ M — v T osillations are also 
discussed. Such an analysis for constraining the mixing angles has been performed in || 
under the approximation of an effective two flavour interpretation of the atmospheric 
neutrino problem either in the — v e or — v T oscillation mode, instead of a full 
three flavour investigation. A detailed analysis combining the accelerator, reactor, solar 
and atmospheric neutrino data had been carried out earlier (pre-LSND) |10| taking a 
different spectrum for Am 2 and assuming the mixing angles to be less than 7r/4. 

The measurement of atmospheric neutrino fluxes is being carried out by the fol- 
lowing groups - Kamiokande, 1MB, Frejus, Nusex and Soudan2. So far, data of most 
statistical significance have been collected by the Kamiokande and the 1MB collabo- 
rations. To reduce the uncertainty in the absolute flux values the usual practice is to 
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present the ratio of ratios R which is defined as, 

D + V»)/(Ve + ^e)obsvd m 
il — -7 _ . ., (1) 

where MC denotes the Monte-Carlo simulated ratio. For neutrinos of energy less than 



~ 1 GeV, 1MB finds R = 0.54 ± 0.05 ±0.12 |11] in agreement with the Kamiokande 
data R = 0.601 

o!o5 i 0-05 in this energy range [§, [12|]. Recently the Kamiokande 
collaboration has published the results of the mesurement of the flux ratio in the 
multi-GeV energy range ||. They found R = 0.57to;o7 ± 0.07 in good agreement 
with the sub-GeV value. All these data show that R is smaller than the expected 
value of unity, a result that might be explained by neutrino flavour oscillation fl"3"| . 
Another aspect of this measurement that can independently point towards neutrino 
oscillation is the dependence of R on the zenith-angle. The multi-GeV Kamiokande 
data reveals a dependence on the zenith-angle unlike the sub-GeV data, though the 
statistical significance of this result has been questioned [fL4 |. For the purpose of this 
paper we use the sub-GeV Kamiokande results. 

The LSND group searches for — > V e oscillations using V e appearance. The V e s 
produce neutrons via the reaction V e p — > e + n which in turn are captured by protons 
producing a 2.2 MeV 7. An excess of beam-on events with a 7 of the above energy in 
time and space coincidence with an electron in the energy range 36 MeV < E e < 60 
MeV is considered as a signal for V e . The mean source-detector distance is 30 metres. 
The initial LSND data reports an excess of 16.4l|^± 3.3 events over the estimated 
background which, if interpreted in terms of neutrino oscillations, corresponds to a 
probability P VfiVe of (0.34^°;?° ± 0.07)%. 

Other appearance experiments searching for — > V e oscillations are KARMEN 
at the ISIS spallation neutron facility |TJJ and the BNL-E776 []nj. These experiments 
are consistent with no neutrino oscillation. KARMEN has so far quoted an upper limit 
on the oscillation probability as Pp %? e < 3.1 x 10~ 3 (90% C.L.) wheareas from the two 
flavour exclusion areas presented by BNL one gets P77 „ e < 1.5 x 10~ 3 (90% C.L.). In 



ref. JI| the LSND group has shown that some of the areas allowed by them in a two 
flavour analysis are disfavoured by KARMEN and BNL-E776. In this paper we confine 
ourselves to the LSND data for constraining the paremeters. 

Reactor experiments searching for neutrino oscillation are GOSGEN, KRASNO- 
YARSK and Bugey. These experiments provide exclusion regions in the Am 2 - sin 2 29 
parameter space by non-observance of neutrino oscillation. The maximum exclusion 
is by Bugey which measures the spectrum of T7 e , coming from the Pressurized Water 
Reactors running at the Bugey nuclear power plant, at 15, 40, and 95 metres using 
neutron detection techniques. The 90% C.L. exclusion contour implies l—P Ve p e < 0.05 



G3- 

The general expression for the probability that an initial v a of energy E gets 
converted to a Vp after travelling a distance L in vacuum is 

ixL 

Pv a vp = $af3 ~ 4 I'nil' hl'njl' !j Sm 2 ( — ) (2) 

Ay 

where Ay = 2A7m(E l/ /MeV)(eV 2 /A ij ), A 4i = rrij 2 — m 2 . The actual forms of the 
various survival and transition probabilities depend on the spectrum of Am 2 assumed 
and the choice of the mixing matrix U relating the flavour eigenstates to the mass 
eigenstates. The most suitable parametrisation of U for the mass spectrum chosen 
by us is U = R13R12R23 where i?y denotes the rotation matrix in the zj-plane. This 
yields: 



U 



I C12C13 S12C13C23 - S13S23 C13S12S23 + S13C23 \ 

— s 12 C12C23 C12S23 

V ~ S13C12 — S13S12C23 — C13S23 — S12S13S23 + c 13 c 23 / 



(3) 



where Cij = cos#y and Sy = sin^y. We have assumed CP-invariance so that U is real. 
The above choice of U has the advantage that #23 does not appear in the expressions 
for the probability for LSND and Bugey. We now turn to the implications of the above 
mixing matrix and the chosen mass ranges on the various probabilities. 
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(i) LSND 

In order to see the impact of three neutrino generations, we first note that for the energy 
and length scales relevant for LSND A23 >> L and the term involving sin 2 (7rL/A23) — > 
0. Further, A13 ~ A i2 and (El) simplifies to 



Pv,v e = 4c^s^ 3 sm 2 (7rL/Ai 2 ) (4) 

(ii) Bugey 

For Bugey, the neutrino energy ranges from 2.8 - 7.8 MeV whereas L is typically 
~ 40 metres. Then A23 >> L, so that sin 2 (7rL/A 2 3) — > 0. On the other hand A12 = 
A 13 << L so that sin 2 (7rL/Ai2) and sin 2 (7rL/Ai3) average out to 1/2. Then the relevant 
probability is 

PjfeUe = 1 — 2c 13 c 12 + 2c 13 c 12 (5) 

(iii) Atmospheric neutrinos 

In a three flavour mixing scheme ([3]) is given in terms of the neutrino transition and 
survival probabilities as 

R= P ^ 

1 VeV e T r . MC 1 i/^Ue 

where tmc = (^e + v e )/(v tl + V^) as obtained from a Monte-Carlo simulation. Notice 
that for neutrinos in the energy range ~ (0.1 - 1) GeV travelling through a distance 
ranging from ~ (10 - 10 4 ) km, A12 = A13 << L and sin 2 (7rL/Ai2) and sin 2 (7rL/Ai3) 
can be replaced by their average value 1/2. Taking this into account, the probabilities 
appearing in (||) can be expressed as 

Pv e v e = l-2c 2 3C 2 2 + 2c 4 3C 4 2 -4(ci3Si 2 C23-Si3S23) 2 (Ci3Si2S23 + Si3C23) 2 sin 2 (7rL/A 2 3) (7) 
P V/tVe = 2c 2 3 C 2 2 S 2 2 - 4c 2 2 C 23 S23(ci3Sl2C23 ~ S13S23) (C13S12S23 + Si 3 C 2 3)sin 2 (7rL/A23) (8) 

p u^ = 1 - 2c? 2 s 2 2 - 4c 4 2 c 23 s 23 sin 2 (7rL/ A 23 ) (9) 

The results of the combined analysis of the above three experiments are presented 
in figs. 1 and 2 in the large s 2 3 and small s 2 3 limits respectively. It is sufficient to 



consider these limits as the allowed values of sf 3 are confined in these ranges. As seen in 
(Hj), the parametrisation chosen for the mixing matrix U ensures that the probability 
relevant for LSND is independent of the mixing angle S23. From the value of Pv^v^ 
quoted by the LSND group |1J one can find the allowed area in the s\ 2 - s\ 3 parameter 
space for fixed values of the ratio Am 2 L/E. The following constraints are found: 
for S12 very small (~ 0) or very large (~ 1), S13 ranges from < S13 < 1 while for 
intermediate values of Si 2 , only very large s± 3 values are allowed. This is between the 
solid lines in fig. 1 (2) for large (small) values of S13, in the limit of sin 2 (7rL/Ai2) ~ 
1. From eq. (|5]) the probability for Bugey is also a function of the same mixing angles 

512 and S13 only, so that, using their result one can further rule out a portion of the 
parameter space - namely, intermediate S13 values at small S12 - which were allowed 
by LSND. This is shown by the dashed lines in figs. 1 and 2 implying the following 
limits for small s 2 l2 {< 0.0018): either s 2 13 >~ 0.97 (fig. 1) or s 2 3 < 0.026 (fig. 2). In 
the other regions of the parameter space the LSND data puts more severe constraints 
than Bugey. 

Our approach next is to determine how much of the combined allowed area from 
LSND and Bugey is consistent with the atmospheric data for fixed values of S23. The 
sub-GeV Kamiokande data implies 

0.48 < R < 0.73 (90%C.L.) (10) 

Imposing this constraint, one finds that the allowed parameter space (shown shaded 
in figs. 1 and 2) depends on the chosen s 23 . In general there are three regions: 

(i) The large sf 3 (>~ 0.97) - small s^ 2 (<~ 0.1) region shown in fig. 1. In this limit it 
is the — v e oscillation that dominates. Considering the limiting case of S12 — > and 

513 — > 1, the relevant probabilities assume the forms: 
p 1 _ O r 2 2 p o„2 2 p ~ 1 _ o r 2 q 2 

*v e v e — J- zc 23 5 23> r Vn u e — z,c 23' b 23' r v ll v lt — 1 z ' c 23' b 23 

From these expressions it is clear that in this limit P v Vr — 0. 

(ii) The large sf 3 and intermediate s 2 2 zone also shown in fig. 1. To understand the 
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transitions in this range we examine the various probabilities in the limit sf 3 — ^ 1. In 
this limit eqs. (0) - (fj|) become 

P ~ I _ 9r 2 s 2 P ~ 2r 2 r 2 s 2 P ~ 1 — ?r 2 s 2 - 2r 4 r 2 s 2 

-T v e i/ e — -L zc 23 6 23' - r ^M I/ e — z ' c 12 c 23 i '23' r v»Mn — ^ L 12*12 z ' c 12 c 23' b 23 

This is the region where the depletion can be due to both the channels simultaneously 
excepting in the special case of s 23 — > when this reduces to solely — v T oscillation. 
From fig. 1 one also notices that irrespective of the choice of S23, large values of s\ 2 
around ~ (0.85-1) are disfavoured by the atmospheric data. In this region v e — u T 
conversion is effective. 

(iii) The small s\ 2 - sj 3 zone - < s\ 2 < 1.8 x 10~ 3 , < s 2 3 < 0.01 - a look at 
the various survival and transition probabilities reveals that this is a region where the 
depletion is mainly due to — v T oscillation. This can be easily seen by considering 
the limiting cases s 12 , s 13 -> 0, when eqs. (0) - (|) give P VpVe ~ 1, P VfI „ e ~ 0, P VflVfl ~ 
1 - 2C23S23. Substituting these in (|10[) one finds 0.162 < s^ 3 < 0.838. There is a 
sharp cut-off as s 23 crosses 0.162 and for practically all intermediate values upto 0.838, 
the whole of the parameter space allowed by LSND and Bugey is consistent with the 
atmospheric neutrino data. Thus in this regime we show the allowed region for only 
one representative s 23 . We have numerically checked that the allowed region is the 
same as the one presented in fig. 2 for all other s 23 in the above range. 

In our analysis we have fixed A12 — A13 at 6eV 2 , where LSND is most sensitive 
and sin 2 (7rL/Ai2) — * 1, maximising the oscillatory term. As discussed in || it remains 
to be seen what best-fit value, consistent with KARMEN and BNL-E776, emerges when 
more data is accumulated. Our results remain unchanged as long as it is permissible 
to use the above limit. 

For the atmospheric neutrino case we approximate the sin 2 (7rL/A23) factor by its 
averaged value 0.5 as is often done in the context of the sub-GeV data [pl| [L9|, |2(J. 
This approximation can be improved by an averaging over the incident neutrino energy 



spectrum, the zenith-angle of the beam as well as the final lepton energy [18, 1C). tmc 



is taken to be 0.45 from a detailed Monte-Carlo simulation including the effects of 



muon polarisation [^T| . 

Finally let us discuss the implications of the parameter space allowed by the 
Kamiokande atmospheric neutrino, LSND and Bugey data for the z/ M — v T oscillation 
search at CHORUS |22[ and NOMAD [23]. These experiments use the beam from 



the CERN SPS with the mean energy ~ 30 GeV and the approximate source-detector 
distance is 0.8 km so that A23 >> L and 

Py.ur = 444^3 sin 2 (7rL/Ai 3 ) (11) 

With the CERN SPS designed to deliver 2.4 xlO 19 protons, CHORUS and NOMAD 
are sensitive to a minimum oscillation probability of 10 -4 . For A12 = A13 in the LSND 
range of ~ QeV 2 , sin 2 (7rL/Ai 3 ) ~ 0.04, whence (|TT|) is P VyMT ~ 0.16c 2 2 s 2 2 s 2 3 . Then for 
the three allowed regions in the s\ 2 — s 2 3 plane one gets: 

(i) In the large s 2 3 , small s\ 2 zone P Vy/Vr can be slightly greater than 10~ 4 being 
marginally within the reach of these experiments. This is the — u e oscillation zone 
for atmospheric neutrinos. 

(ii) For large s 2 3 and intermediate values of s 2 2 , P v Vt is ~ 10~ 2 , which is well within 
the reach of CHORUS and NOMAD. Recall that this is the genuine three generation 
oscillation regime for atmospheric neutrinos where both — v T and — v e modes are 
operative, excepting for the special case of S23 — for which it is just — u T . 

(iii) In the limit of both s 2 2 , s 2 3 small, P v Vt is very small and this regime, where the 
atmospheric anomaly is due to — v T oscillation, cannot be probed by CHORUS and 
NOMAD. 

For the chosen values of the mass- differences a simultaneous solution to the solar 
neutrino problem is unobtainable unless one invokes a sterile neutrino. Work is in 
progress in this direction ||24|| . 

In conclusion, we have obtained the mixing angles compatible with atmospheric, 
LSND and reactor experiments (in particular Bugey) in the context of three flavour 
mixing. Keeping Am 2 fixed at the best fit values obtained from two generation analyses 
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of the LSND and atmospheric data, we find three regions of parameter space that can 
account for all three experiments simultaneously. Our results differ from an analysis 
presented in || in that we find the mixing angles to be less restricted. Our method, 
which takes into account the possibility that the depletion of the atmospheric neutrinos 
can be simultaneously due to — v e and — v T oscillations, is more general and 
includes the constraints obtained in as a special case. A direct comparison of the 
values obtained for the mixing angles is, however, not proper because the definitions 
of the mixing matrices are different. The sensitivity of the accelerator based neutrino 
oscillation experiments at CERN, CHORUS and NOMAD, is different in the three 
allowed zones and thus they can distinguish between these sectors of the parameter 
space. We find that CHORUS and NOMAD are most sensitive to that part of the 
combined allowed area where the atmospheric neutrino anomaly is due to — v T and 
Vn — u e oscillation modes simultaneously. 
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while A.R. has been supported in part by the Department of Science and Technology 
and Council of Scientific and Industrial Research, India. We thank S. Mohanty for 
sending us some of the useful literature. 
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FIGURE CAPTIONS 



Figure 1: The 90 % C.L. allowed region in the s\ 2 — s\ z plane from LSND is between 
the solid lines, that from Bugey is above the dashed line while the combined allowed 
area including the Kamiokande sub-GeV data is shown shaded. 

Figure 2: Same as in fig. 1 excepting the region allowed from Bugey is below the 
dashed line. 
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